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Abstract —Vehicle-sharing systems are becoming important 
for urban transportation. In these systems, users arrive at a 
station, pick up a vehicle, use it for a while and then return 
it to another station of their choice. Depending on the type 
of system, there might be a possibility to book vehicles before 
picking-up and/or a parking space at the chosen arrival station. 
Each station has a finite capacity and cannot host more vehicles 
and reserved parking spaces than its capacity. We propose a 
stochastic model for an homogeneous car-sharing system with 
possibility to reserve a parking space at the arrival station when 
picking-up a car. We compute the performance of the system and 
the optimal fleet size according to a specific metric. It differs from 
a similar model for bike-sharing systems because of reservation 
that induces complexity, especially when traffic increases. 

I. Introduction 

A. Context 

Over the past decade, vehicle-sharing systems have ap¬ 
peared as a new answer to mobility challenges, like reducing 
congestion, pollution or time of travel. As a consequence, the 
number of cities equipped with bike-sharing sytems is almost 
400 (with 60,000 to 120,000 trips a day in Paris for instance) 
and car-sharing systems are beginning to spread as well. 

Those systems can be described as follows. For a typical 
bike-sharing system, users simply arrive at a station, pick a 
bike if there is one, otherwise leave the system, use it to 
go to another station and leave it there, if there are spaces 
available. If no space is available, the users have to find 
a neighbouring station to return their bike. For car-sharing 
systems, users have the possibility to reserve either a parking 
space only or both a car and a parking space in the destination 
station before picking-up the car. Again, they can only do this 
if there are cars and spaces available. Otherwise they leave 
the system. This basic description could be refined, especially 
users, instead of leaving the system unhappy, could look for 
a neighbouring station with available resources. The lack of 
resources is one of the major issues which oblige operators 
to perform regulation to maintain the reliability of the service 
against other transportation modes (see H2I and ffl). 

B. Related works 

Due to the success of vehicle sharing programs and the 
crucial problem of allocation of ressources they face, both 
vehicles and parking spaces, there have been many studies on 
redistribution in bike-sharing systems. See |2J and reference 
therein. Much of the work concerns optimization applied to 
deterministic models. 


Stochastic approaches have been developed recently for 
bike-sharing systems. For pioneer papers, see II and II 
for a model with stations with infinite capacities, where the 
problem of lack of parking spaces is avoided. In 0], the first 
homogeneous model addressing the problem of the availability 
of parking slots aims to investigate the influence of parameters 
such as demand or capacity on the performance of the system, 
defined in terms of the proportion of problematic stations 
(empty or full stations). The paper 0j deals with a set of 
N identical stations with capacity K. Users arrive at rate A 
in each station, pick-up a bike if there is one, join the pool of 
riding users and then return their bike in a station chosen at 
random, after a trip duration with exponential distribution with 
parameter //. If the station is saturated, they reattempt with the 
same strategy until they succeed. The authors prove that, in 
this model, the minimum stationary proportion of problematic 
stations as the system size gets large is 2 j(K + 1). This 
minimum is reached when the average number of bikes per 
station s equals K/2 + A /p. 

The paper 01 also assesses the consequences of different 
incentive policies such as choosing when returning the bikes 
between two stations the one with less bikes, which drastically 
reduces the proportion of problematic stations. In another 
paper 0, they show how to extend the results on the basic 
model to a inhomogeneous framework, considering subsets of 
stations with the same parameters. In those papers, the Markov 
process considered is the empirical distribution of the number 
of bikes in stations, i.e. a vector with the proportion of stations 
with k bikes at time t. The mean-field limit gives an ODE 
satisfied by this empirical distribution when the number of 
stations N tends to infinity. Then, it is proved that the solutions 
of this ODE converge with time to a unique equilibrium point, 
and that this equilibrium point is the concentration point of the 
invariant measure of the empirical distribution process as the 
system gets large. It gives here that, in both homogeneous and 
inhomogeneous cases, the number of bikes in a station follows 
a truncated geometric distribution at equilibrium, when the 
system gets large. 

In an upcoming paper |6), Tibi and Fricker prove directly 
the convergence of finite marginals of the invariant measure of 
the Markov state process, i.e. the vector of length N containing 
the number of bikes of each station at time t. Thanks to 
irreducibility and reversibility, it has a product form and, using 
a local limit theorem, this gives moreover that the number of 
bikes at a fixed number of stations are independent in the limit 
N —> oo. This property is always difficult to reach via the 


mean field method. An interesting result, original as far as we 
know, in J5) is that the Markov process including both stations 
and routes states has still a product-form invariant measure, for 
the model with fixed size. It extends the result of George and 
Xia@, for the scaling considered along our paper, where both 
numbers of stations and bikes tend to infinity, and for finite 
capacities. 

Waserhole et al. lfl2ll propose a Markovian model for car¬ 
sharing systems where they use a fluid limit approach to be 
able to maximize the average profit of the operator through 
revenue management techniques. The paper ED investigates 
the impact of reservation through different policies in vehicule- 
sharing systems. 

C. Outline of the paper 

To our knowledge, our paper presents the first stochastic 
analysis of a large-scale car sharing system. 

We present an homogeneous model with M cars and N 
stations with capacity I\ where users have to reserve a parking 
space in a destination station chosen at random when they 
pick-up a car. Users arrive at rate A in each station. If they 
find a car and if there is some free space in their destination 
station, they pick-up a car and make the reservation at the 
same time. Otherwise, they leave the system. 

Unlike bike-sharing models, neither the empirical mea¬ 
sure process nor the state process are reversible. Even though 
we know that there exists an invariant measure for a fixed 
N for these irreducible finite-state space Markov processes, 
it is quite untractable. Thus the aim is to obtain large-scale 
asymptotics to understand the large-scale behavior of the 
system. It means that M and N are large with the number of 
cars per station tending to a constant s. This sizing parameter 
is a key parameter of the system. For that, we are still 
able to adapt the mean-field limit approach in 13 when the 
system gets large. The main difference is that, due to the 
reservation, a station is described by two components: the 
vehicules and the reserved parking spaces which makes the 
analysis of the equilibrium more tedious. Indeed, recently, 
large-system asymptotics have been successfully applied in 
many contexts in com- munication systems (see 0, |8|] and 
others). But it is, as far as we know, the first system where 
the underlying process gouverning the dynamical system is 
not one-dimensional. 

The analysis of the equilibrium point gives the following: 
when the system gets large, the steady-state numbers of 
vehicles and reserved spaces at each station have a product 
form distribution with a geometric and a Poisson terms on a 
constraint space. The two parameters involved are solutions 
of two fixed point equations. This allows us to study the 
behavior of the system. As intuition suggests it, we prove 
that all stations are problematic when the average number of 
vehicles per station equals the capacity, which makes quite 
a difference with the homogeneous model for bike-sharing 
systems, or when there are no vehicles. We investigate the 
fleet sizing problem. We obtain asymptotics in the two cases 
of light and heavy traffic. We prove that in light traffic case. 


reservation has little impact on performance, unlike the heavy 
traffic case. The main difference with bike-sharing systems is 
that the best performance degrades with traffic and that the 
corresponding fleet size remains under the overall capacity. 

Then, a model with reservation of both resources, vehicles 
and parking spaces, is proposed. An approximated simple 
model is studied and our simulations show its relevancy. The 
analysis shows that the performance of such a system behaves 
as in a system with simple reservation but with more traffic. 

Section m describes the model, the Markovian process 
and the limiting ODE. In Section [Till we prove the uniqueness 
of the equilibrium point of the ODE. In Section lfVl we analyze 
the equilibrium and give performance results. Section M deals 
with the reservation of both resources. Then Section eh 


concludes the 

paper. 

II. Stochastic Model 

A. Main Notations 

N 

Number of stations. 

M n 

Total number of vehicles. 

Sn 

Average number Mn/N of cars per sta¬ 
tion. . 

K 

Number of parking spaces in a station, 
also called capacity of the station. 

A 

Arrival rate of users at a station. 

1/p 

Average trip time. 

v?(t) 

Number of vehicules in station i at time 

t. 

Number of parking spaces reserved in 
station i by users travelling at time t. 

R?(t) 

y k N i(t ) 

Proportion of stations with k cars and l 
reserved parking places at time t. 


Limit of Y*(t) as N tends to infinity 
(described by an ODE). 

y 

Equilibrium point of the corresponding 
ODE. 

X 

Space of station states. 


B. Model description 

The system is a set of N stations with capacity K with 
Mjv vehicles. As N tends to infinity, sn = Mn/N tends 
to s called the average number of vehicles per station. At a 
given station with capacity K, users are supposed to arrive 
with rate A. An arriving user at station say i has a destination 
say j chosen at random. If there is no car available in the 
station origin i or if there is no available parking space in 
the station destination j , the unhappy user leaves the system. 
Otherwise, she picks-up a car in station i and simultaneously 
makes a reservation in station j. Then, the journey between 
station i and station j takes an exponentially distributed time 
with mean 1 /p. The user returns her car at station j and leaves 
the system. 

Note that unlike bike-sharing systems (S3)- users do not 
have to look for a station to return their vehicle until they find 
a parking space available, thanks to reservation. 


In this paper, we focus on this homogeneous model. But 
we are able to easily extend the result to a heterogeneous 
model consisting in a finite number of clusters with for a 
station of cluster i a capacity K,, an arrival rate A ? , and 
probability pi of choosing a destination in cluster i, instead of 
1/N. 

C. The Empirical Measure Process 

Let us define x = {(A:,Z) £ N 2 ,fc + l < K} and let 
Y k j (t) be the proportion of the N stations with k vehicles 
and l reserved parking spaces at time t, that is 




1 

N 


N 

' 52 1 {V i N (t)=h,R?(t)=l} 

2=1 


where (t) is the number of vehicles and it) the number 
of reserved parking spaces at station i at time t. If 'P(x) is 
the set of probability measures on y, let 

y = {y €V{x), ^2 ( k + l )Vk,i=s} 

(k,l)£x 

and 

N 

y N = {y £ V{x),yk,i e Jj, ( k + l )yk,i = sn}- 

( k,l)Ex 

As the system is homogeneous, the process ( Y N (t )) = 
i ex is a Markov process on y N , described as 
follows. Suppose the process (Y N (t)) is at state {yk,l)(k,i)ex- 
There are two different types of transitions: 

• Cars picked up. The arrival rate of users in a station in 
state (fc, l) is A y k jN if k > 0. A user, who arrives at a 
station, makes at the same time a reservation in a station 
in state (k',l') with probability yk',v if k! + V < K. 
Therefore, the transition rate is of A y k 'py k yN if k > 
0, k' + l' < K. The arrival causes y k j to decrease by 
1/N and y k -iy to increase by 1/JV. And the reservation 
causes y k > p to decrease by 1/N and yw ,v+i to increase 
by 1/N. 

• Cars returned. When a car arrives at its reserved parking 
space in a station in state (k,l), y k ,i decreases by 1/At 
and yk+i.i-i increases by 1/At, as the reserved space 
is replaced by a car. Considering stations in state (fc, Z), 
the number of reserved parking spaces is lNy k y. As the 
trips are exponentially distributed with mean 1/p, this 
transition occurs at rate plNy k y. 

The jump matrix Q N of the process ( Y N (t )) is thus given 
by, for y, y' £ y N , 

f Q N (y, y') = a yk’,i>yk,i if fe > 0, k' + T < K 

[Q N (y,y+ jf( e k+i,i-i - e fc ,z)) = plNy k ,i 

where y' = y + ^-(e fc _ u - e k y + e k p V+1 - e k p V ) and 
i e k,l)(k,i)ex are the vectors of the canonical basis of RM. 

It is thus irreducible and on a finite set y N , which implies 
it admits an invariant measure denoted by fl v . As I l Y is 
analytically untractable, we use a mean-field limit when N 


tends to +oo to give a large-scale asymptotic. For that, we 
show that (Y N (t))o<t<T tends to y(t)o<t<T solution of an 
ODE. Then, in Section [Till we prove that the ODE admits a 
unique equilibrium point y. This equilibrium will be analyzed 
in Section Hvl 

D. Dynamical System 

By standard arguments, for t > o, (y JV (f)) te[0)T] 
converges in distribution to (y(f))te[o,T] unique solution with 
y(0) fixed of the following ODE 

y{t) = yk,i(t) \ yk,i{t)pl(ek+i,i-i-e k ,i) 

( k,l)ex \(fc',i') 6x 

+^yk',l'lk>0;k'+l'<K(ek-l,l — e fc,Z + e fc',I'+l — . (1) 

The first term corresponds to the rate at which users 
return cars at the reserved parking space, and the second term 
corresponds to the rate of simultaneous arrival and reservation. 
Let us introduce some additional notations. Let 

K-k 

ys(t)= M*). w,(t) = E ( 2 ) 

k+l=K 1 =0 

be respectively the limiting proportion of saturated stations at 
time t and the limiting proportion of stations with k vehicles 
at time t (0 < k < K). Then, by splitting the second term of 
its right-hand side, the ODE (|TJ rewrites 

y{t)= ^2 yk,l(t) (A(l — ys)i^k-i,l ~ efc,;)lfe>o 

(k,l)ex 

+A(1 — t/o,.)(e/c,;+i — & k ,i)1- k +i<K + pK^k+iy-i — £k,i)) ■ 

of the form 


y(t) = y{t)L v (t) (3) 

where y{t)L y ^ is the product of the vector y{t) by the jump 
matrix L y (t) defined in the following section. This means that, 
when N tends to infinity, the empirical distribution y(t) of 
the stations evolves in time as the distribution of some non- 
homogeneous Markov process on y, whose jumps are given 
by L y (t), updated by the current distribution y(t). 

III. Steady-State Behavior 
A. Probabilistic Interpretation of the ODE 

The jumps of the Markov process with generator L y on 
X = {( k,l ) £ N 2 ,k + l < K} defined in Section Hlj can 
be seen as those of the number of customers in two coupled 
queues (see Figure Q}: a typical station can be described as 
a tandem of two queues: a M/M /oo queue for reservations, 
with arrival rate A(1 — 2/o,.) an d service rate p, and a M/M/1 
queue for vehicles, where customers come from the former 
queue, with service rate A(1 — ys)- Moreover, this system is 
a loss system: the total number of customers is less than or 
equal to K. Then the problem is now to find the stationary 
measure of this system of two queues in tandem. 


R: M/M/ oo 


A(1 - 2/o„) 


P 


V: M/M/l ^ (i- K) 
-> 


Fig. 1. A typical station as a tandem of two queues with overall capacity K. 


Let pn be the arrival-to-service rate ratio of the M/M/oo 
queue, and pv of the M/M/1 queue, defined as 


pr(v) 

Pv(y) 


-(! - 2/o,.), 

P 

1 - 2/o,. 

1 - ys ' 


(4) 

(5) 


In this case, it is well-known (see HD for example) that 
the invariant probability measure n(y) associated to L y has a 
product form, given for (k,l) £ X by 


nk,i{pv(y),PR{y )) 


1 

z(pv(y), pr(v)) 


pn/y) 1 

i\ 


Pv(y) k 


where Z(p v (y), Pr{v)) is such that 


(6) 


E n k ,i(pv(y),pR(y)) = i- 

(k,i) ex 


With a slight abuse of notations, we will use both n (y) 
and 7r (pv{y), PR(y)) in the following. 

Because n(y) £ y, we also have 

s = ^2 ( k + l ) 7r k,i(y)- (7) 

(fc.Oex 


which is equivalent to s = E( V+R) where (R, V) is a random 
variable with distribution 7r(y). 

To investigate the steady-state behaviour of the model, 
we study the equilibrium points y of the ODE 0. Such 
equilibrium points y are the y satisfying 


v{pv(y),PR{y)) =y ( 8 ) 

and equations 0. (0 and 0. Therefore, it is also equivalent to 
finding couples ( pv,PR ) such that y = n( pv,Pr) satisfying 


PR 

= -(1-7 t 0j Xpv,Pr)) 

P 

(9) 

PV 

1 — k 0 ,Xpv,Pr) 

1 - tts(pv,Pr) 

(10) 

s 

= E {R + V ) 

(11) 


where (R, V) is a random variable with distribution 

k(pv,Pr)- 

In conclusion, the question of finding a measure on \ is 
reduced to find a couple of real numbers ( pv,PR )■ The two 
following sections prove the existence and uniqueness of the 
equilibrium point. 


B. Existence of an Equilibrium Point 

Due to the Section ITlI-AI we know that finding an equilib¬ 
rium point y of the ODE 0 is equivalent to finding y solution 
of equation 0. Since s < K, for each y £ y, ys < 1. 
Therefore, by definition (see equations 0, 0 and 0), 
V l— ^ (Pv{y): Pr( y)) is continuous from the compact and 
convex set y to itself. Thus, the existence of y is given by 
Brouwer’s theorem. 


C. Uniqueness of the Equilibrium Point 

Compared to the bike-sharing model in 0, reservation 
makes calculations much more tedious. But we still prove the 
uniqueness of the equilibrium point. Let us present the main 
steps of the proof. 

First, we prove that Equation m is true for any 
( Pv,Pr )■ Second, we consider equation 0 as an implicit 
function between pr and pv which gives pv as a strictly 
increasing function of pr. The third step consists in studying 
the monotony of the right-hand side of equation CD as 
a function of pv and pr (as for bike-sharing systems, it 
increases), and then, in using the diffeomorphism between pr 
and py to conclude on the uniqueness of the solution of 0 
and (ITTb . 

From now, we assume © to be satisfied since we 
know there exists solutions. We need to find the solutions 
( Pv,Pr ) °f equations 0, (ITOl) and dTTT) . Let us first notice 
that necessarily ( pv,pR ) £ T = [0,oo[ x [0, X/p[ according 
to equations © and (QB. 

The partition function Z(pv,Pr) of the invariant proba¬ 
bility measure plays an important part in the following, so we 
mention some of its properties. 

Property 1: (i) For any ( Pv,Pr ) £ T = [0, oo[x [0, \/p[. 


Z{pv,Pr)= ^2 

( k,i)e x 


Z 


dz 

dpR 


(1 -7T S )Z, 


( — ,Pr ) Pv = Z(pv, PrPv) if Pv > 0. 
\Pv ) 


(ii) For any ( pv,PR ) solution of equation 0, 


(/i ~ PR ) Z ( pR,Pv ^ 


A 

P 


K 


E 

1=0 


2r 
l \' 


Proof: Simple algebra gives the result. For the third 
equation of (z), using the change of indexes k! = K — k — l, 



4= E = 

(k,l)& X i= 0 fc=0 

K i K—l i 

E Pr M+i __ 'X' Pr n k'+i 

~/y Pv - 2 ^ ~jpPv 

l— 0 k'= 0 (k',l)£X 

= Z(p Vl p R pv). 























Property 2: Equation (flOb is true for all (py, Pr ) S r = 
[0, oo[ x [0, X/p[. 

Proof: Using definitions given by (0, 


Pv( 1 - k s (pv,Pr)) = 


1 


z (pv,Pr) 


Pv 


k+l<K 


E 


1 


E 


Pr D k 

Z(p V ,p R ) ^ l'. v 

\yv ) yni k>o k+ i^ K 


= 1 — n 0 ,.{pv, Pr)- 


Thus, the second term of the right-hand side of equation ( fi~3l ) 
can be rewritten 



PvPr k 

(K — k)\ 


\ 1 K J K n k n K-k 

A 1 Pr ST' PvPr 

P Z{Pv,Pr) " P- 

\ K K i—j j 

A Pr ST' Pr Pv 

pZ(p v ,p R ) i\{K — j)\ 


(15) 


Property [2] ends the first step of the proof: It allows us 
to focus on equations 0 and (fill only. It means that the 
set of equations 0, (flOt and (fill is equivalent to the set of 
equations 0 and (fill . 

Theorem 1: There exists a strictly increasing diffeomor- 
phism cj> : [0, A /p[ —> [0, oo[, and ip = such that (py, Pr) 
is a solution of equation 0 if and only if py = <t>(p R ). 


A pf f V fcfPy V \ 

PZ{PV,PR) \o <kt<K i[ ( K ~ j)l h iKK - j V) 

For the first term in the sum of the right-hand side of 
equation (IT5l ). 


E 


o <j<i<K 


l—J 7 

Pr Pv 

i\{K - j)\ 


1 

K! 


E 

(k,l)£x 


K'-P l R P k y 


(k + l)\{K-k)\ 


> 


Z(pv, Pr) 

I<\ 


Proof: It uses the implicit function theorem. Due to 
Property |T](m), (py, p R ) G 1 is solution of equation 0 if and 
only if f (py , Pr) = 0 where / is the C°° function defined by 

Kpv,Pr)=( - Pr]z(p r ,p v ) - Elf ' ( 12 ) 

\P J P l=0 <• 


using that, for any (k,l) £ 

(k + l)\(K-k)\ ~ n k + l-i - Pf 

Moreover, the second term in the sum on the right-hand side 
of (fl5T > is positive. Thus equation (IT5l) yields that 


And we already know by Section UlI-BI that there exists 
a solution of equation 0. Let ( pv,p R ) be such a solution. 
To prove the existence of a strictly increasing diffeomorphism, 
we have to show that 




First, using Property [T| ( i ), 


df_ 

dpR 


(pv,Pr) = -Z(p v ,pR) 


A 

- PR 

p 


E 

k+l<K -1 


A 


K-l 


E 

/—0 


Pr 

l ! 


(13) 


thus, substracting equation l lT3l > to equation ( IT2l >. 


f(Pv,pR.) 


-^-(pv,Pr) = Z{p v ,p R ) 

opR 



n k n K ~ k 

PvPr. 

(. K-k)\ 



From equation (fl2l >. while f(py , pr) = 0, it holds that 


\ \ K i 

-~ PR = _-_ yE*. 

P P z {pr,Pv) “ P 



PvPr k 

(. K-k)\ 


Apf 

/i K \' 


> 


(16) 


Additionally, using that Z > 0, we can deduce from 
equation (fl4l > that f(pv,PR) — ■§^{pv,Pr) > 0. Therefore, 
while f(pv,PR.) = 0, 

df 

-R ( PV,PR) < 0. 

OPR 

It is easy to see that ~§~{pvi Pr) is positive: Since p R < 
A /p and Z strictly increasing in py , 

df , , ( A \ dZ . . „ 

d^ {pV ' pR) = [p ~ PR ) d^ {pV ' pR) > °- 

One can easily check that if py = 0 then p R = 0. 
Similarly, when py tends to infinity, p R tends to A/// to keep 
/( Pv,Pr) = 0. It ends the proof. ■ 

Theorem Q] means that the solutions of equation 0 can 
be expressed with only one parameter. It will also be useful 
for a further study of equilibrium as it allows to differentiate 
expressions. This concludes the second step of the proof. 

Once we have used equations 0 and (fTot . let us focus 
now on equation (fTTV We prove that E(/f + V) as a function 
of py and p R is strictly increasing in each parameter, and 
then using Theorem Q] we conclude to the uniqueness of the 
equilibrium point. 

Property 3: The average number s(py, p R ) = E(i? +V) 
of vehicles and reserved places per station, where (R, V) is 


























Fig. 2. Diffeomorphism between py and the average number of vehicles s 
with K = 4. s is strictly increasing with traffic even though the stations’ 
capacity is finite. 


a random variable with distribution tt(pv,Pr), is a strictly 
increasing function of both py and p R . 

Proof: See Appendix IVIII ■ 

The following result concludes the proof of the unique¬ 
ness of the equilibrium point for the ODE. 

Theorem 2: For any s > 0, there exists a unique 
( pv,pr ) solution of both equations (O and (ITTl i. 

Proof: We have to prove that the function from [0, +oo[ 
to [0 ,K[ which maps py to sx(pv,^{pv)) is a strictly 
increasing diffeomorphism. It is true, using both Theorem |T| 
and Property [3| and that 


S K ~ 


ds K ds K 
dpy dp R 




It ends the proof. ■ 

Figure [2] represents the diffeomorphism between s and py 
for different values of traffic A///. As we could have expected, 
it is globally increasing with traffic. 

The proof of the uniqueness of the equilibrium point of 
the ODE is much more tedious than in the model for bike¬ 
sharing systems in E) since we have two parameters related 
by just an implicit function. 

This equilibrium point also allows us to study the per¬ 
formance of the system in Section [TV] We can parametrize all 
quantities with one parameter only and also differentiate them. 


D. Convergence to the Equilibrium Point 

Both the convergence of the dynamical system to its 
unique equilibrium point and the convergence to the steady- 
state empirical measure process to the same point are out of 
the scope of the paper. These questions are still open. In the 
bike-sharing model of S). thanks to reversibility, a Fyapunov 
function is found. Here, the stochastic model for car-sharing 
systems is much more complex: The state process of a tandem 
of queues is not reversible. Recall also that Tibi 0 proves the 


convergence of the invariant measure of the state process by 
its explicit product-form. It uses a local limit theorem. For the 
car-sharing model, no closed-form expression of the invariant 
measure can be expected. 

We thus study y as the heuristic limit of both the dynam¬ 
ical system and the stationary empirical measure Y N ( oo). It 
gives the limiting stationary behavior of our model. 

IV. System Performance 

To manage the system, it would be interesting to know 
how to choose the key parameters K and s to satisfy most 
users of the system, or to maximize profit. In this section, 
we discuss the choice of a metric to assess the homogeneous 
system performance and then give results or conjectures about 
optimal behavior deriving from the main metrics. 


A. Choice of the metric 

By Section [HI] admitting the convergence as the system 
size gets large of the invariant measures, the limiting stationary 
distribution of the joint numbers of vehicles and reserved 
places in a station is y. We mainly focus in our paper on 
a quality of service indicator, called the limiting stationary 
proportion of problematic stations. 

Definition 1 (Problematic Stations): Fet y be the unique 
equilibrium point of ODE ©. The stations with either no 
car or no parking space available are called problematic. The 
limiting stationary proportion I), of problematic stations is 
given by 


Pb — — 0 or I + R — K) — 2/o,. + Us ~ Vo,k 

where (V,R) is a random variable with distribution y. For 
simplification purposes, we also define P b + by 

p b =v o„ +Vs- 

Note that Iff > Pb. We can also expect that Pjf is a 
good approximation of Pb- This metric gives a convenient 
way to compare the performance with the bike-sharing system 
performance described in ED- It is also close to the proportion 
of unsatisfied users. Indeed, a user is satisfied if she can 
pick up a car at her chosen source station and can return it 
at her chosen destination station. In a homogeneous system, 
and if stations are chosen uniformly at random, this occurs 
with probability (1 — y 0 )(1 — y s ). Hence, the probability of 
unsatisfaction is U = y 0 + Vs ~ Vo Vs- Generally, the third 
term is small compared to the others. In our case, we can note 
that Pb > U because of equation (1 1 6k key argument in the 
proof of Theorem [I] which rewrites, using equation (IT5l) . 


K 


Z(pv,Pr )“ 


\ ' Pr 

l\ 


K n k n K ~ k n K 

E PvPr > Pr_ 
k _ 0 (K-k)\ - K\ 


which gives, by dividing by Z(p v ,p R ), y 0 y s > y 0 K . 

The following equation, straightforward consequence of 
equations © and (TITTb . will be useful. 


D + _ o PPR PPR 
b A A py 


(17) 











Therefore, the limiting proportion of problematic stations 
Pb is close (and an upper bound) to the limiting proportion of 
unsatisfied users who can not find a car or return their car in 
the desired stations, which is a frequently used performance 
metric, and P b + is a useful upper bound for both of them. It 
is summarized in the following proposition. 

Proposition 1 (Metric comparison): 

n + >p b >u. 

One might also be interested in the number of trips 
achieved per unit of time. Note that the operator revenue is 
more or less proportional to this quantity. 

Definition 2 (Number of Successful Trips): Let y be the 
unique equilibrium point of ODE ((3}. The limiting stationary 
number of successful trips per unit of time is denoted by T 
and given by 

T = A(1 — U) 

and moreover, T > A(1 — Pb) and approximated by A(1 — P&). 

Note that in this model, the average sojourn time does 
not depend on the state of the stations since the parking space 
is reserved at the destination. Thus, we will not consider any 
metric related to this quantity unlike what is done for bike¬ 
sharing systems in (4). 

In the following, a value of parameter X referred as 
optimal must be understood in the sense of minimizing the 
proportion of problematic stations for a given traffic. Such a 
value is denoted by X*. 

B. Optimal Fleet Size 

We study in this section the influence of parameters on 
performance, especially the fleet size given by parameter s. 

Due to the results of the previous section, mainly The¬ 
orem HI giving the existence of t/>, y function of (pv,Pr), 
and other derived functions as performance metrics, can be 
considered as a functions of py only, because y is solution of 
equation ©. With abuse of notation, these variables will be 
denoted by y(py ), Pb(pv) and s(py). Thus the proportion of 
problematic stations, as a function of the fleet size, is given 
by the parametric curve 

pv ( s(p v ),Pb{pv ))■ 

This curve will play an important role in the following, giving 
the behavior of the system according to the fleet size parame¬ 
ter. It is numerically plotted in Figure [3] We emphasized that 
this curve is not obtained by simulations. The performance 
metric chosen here is compared to the same for bike sharing 
systems (see El for details). 

Further results will rely on an interesting property of 
symmetry as follows. 

Property 4 (Symmetry): For all py > 0, 

4>(1/Pv) = — tp(pv), 

Pv 

Vs^/pv) =Vo Jpv), 

P b (l/pv) = Pb(pv) and P+(l/p v ) = P^(pv). (18) 


Proof: Let py > 0 and pr = ip(py). We also define 
PR = i’i 1 /Pv)- Then, 


-n, ) 1 1 V- (PvPn) K ~ k 

Vs Pv Z(l/pv,p R ) p* E (K — k)\ 

But, with Property Q] Z(l/p Vl p R )p$ = Z(p v ,p v p R ). 
Hence, 


VsiX/pv) 


1 

Z(.Pv,PvPr) 


K 


E 

1=0 


(pvpr) 1 

l! 


Using TheoremQ] since (py, PvPr) is solution of Equa¬ 
tion ©, it holds that pypR = pR , which is equivalent to 
V>(1 /pv) = 1 /pyip(pv). Thus y s (l/pv) = y K {p v ). 

Then, with the results above, 




Pr(pv) k 

PvZ(^,Pr(^))K\ 


= yo,.(pv) + ys(pv) 


pr(pv) k 

Z(pv, Pr)K\ 


Thus, Pb(l/py) = Pb(pv). The same is true for P b + . ■ 

Figure [3 seems to illustrate the presence of a minimum 
for Pb as a function of s. But this result is unreachable and 
conjectured here. Nevertheless, we propose the following weak 
result. 

Proposition 2: For any A//./ > 0, I), has an extremum 
P b * for p v = 1. 

Proof: Existence is clear since Pb is continuous and, by 
definition, P&(0) = 1 and Pb tends to 1 as py tends to +oo 
using equation d 

The value py = 1 comes also from the symmetry proved 
in Property [4] Indeed, using equation (fl8l) . 


Pb(Pv) = - -T P b( l /Pv)- 

Pv 

Therefore, Iff 1) = — P^( 1) = 0 and the sign of the derivative 
changes at py = 1. Hence, for py = la local extremum is 
reached for Pb , and the same holds for P b + . ■ 

Uniqueness of the extremum should come from the 
convexity of the parametric curve py i—>■ (s(py), Pb(py)). 
This yields that it is a minimum. Nevertheless this convexity 
result is still to be proved because the implicit relation between 
py and pr (see Section ED makes calculations tedious. 

Conjecture 1: For any X/p > 0, Pb has a unique mini¬ 
mum Pb* reached for py = 1. 

Conjecture [T] is similar to the result for bike-sharing 
systems (see El for details), that the minimum is reached in 
py = 1. In Figure [3 we observe the existence and uniqueness 
of the minimum conjectured in Q] The curves plotted for both 
systems are convex with a unique minimum. 

Moreover, for the influence of the fleet size on the system 
behavior, the fact that Pb tends to 1 as py tends to +oc makes 
an important difference with the homogeneous bike-sharing 
system. Without reservation, Pb goes to 1 when s tends to 
infinity whereas reservation implies full saturation when s = 








K. This fact will be developed in the following with the study 
of the heavy-traffic case. 

Property 5: For all py > 0, the approximated proportion 
of problematic stations Pp(pv) is increasing with traffic load 
X/p. 

Proof: Define a = A/p > 0 and f a = f/a where / is 
given by equation (IT2t . It holds that 

K f 

fa(pv,PR) = (1 - Pr/o)Z(p V ,Pr) - ^2 (19) 

1=0 

Denote by ip a the diffeomorphism ?/; defined in Theorem Q] 
Let a > 1 be fixed. First let us prove that 

ail) a > Ip cm- (20) 

It is equivalent to prove that faa(pv,apii) < 0 where 
(pv,Pr) different to (0,0) satisfies f a (pv,PR.) = 0, i.e. 
ipa(pv) = Pr- For that, using equation ( IT9t . it holds that 


Z{pv, PR)faa(pv,apn) = 

l K l 

Z(pv,PR)(l~PR/a) ^2 Pv^-a l - z (Pv,PR)^2^-a l . 

( k,l)E X ' 1=0 L ' 

K , 

Since f a (pv,PR) = 0, Z(p v ,p R )( 1 - p R /a) = £ if • 

1=0 

Hence, 


Z ( Pv , Pr) fa.a (pv i &Pr) 


K 


K 


= £? £ Ef* 1 

1=0 (k,l)£x (fe,l)Gx 

K K-kK-k j+l 
= ££ %-(“■-“ j ) 


/c—0 /—0 j=0 


j+Z 

EE £ -o't 


fe=0 i=0 j=JC-fe+l 


( 21 ) 


The hrst term of the right-hand side of the equation (|2TT > equals 
zero as it is symmetrical in l and j. And because j > l and 
a > 1 the second term of the right-hand side of equation i lTil 
is negative, which concludes the first part of this proof. 

Then, by equation (fTTI ). 


jt>+ / \ 0 1paa{pv) 1paa{pv) 

Pf{aa) = 2-. 

aa aapv 

Plugging equation (l20l > in it, we can conclude that P b + (aa) > 
Ppp(a). It ends the proof. ■ 

Property 0 is not so interesting in practice since the 
operator does not have directly access to the parameter pv- 
However, it could be a first step for proving the following 
Conjecture [2] 

Conjecture 2: For any s £ [0. K\, the proportion of 
problematic stations P fe + (s) is increasing with traffic A/p,. 

We conjecture both Property [3 and Conjecture [3 are valid 
for Pi, as well. In fact, assertions given in Conjectures [2] and Q] 


for Pf, can be observed on Figure [3 where I), is plotted as a 
function of s. 

Conj ectu re [3 wou I d imply the behavior of the system with 
reservation is very different from systems without reservation. 
Indeed, for bike-sharing systems, increasing traffic does not 
systematically worsens the situation (e.g. if s is close to K as 
we can see on Figure [3 (see also ID). But in both systems, 
increasing the capacity improves the system performance. 


(a) Capacity of stations: K = 5 



Number of vehicles per station : s 
(b) Capacity of stations K = 10 



Number of vehicles per station s 


Fig. 3. Parametric curve py > (s,P&) for K = 5 (a) and K = 10 (b). 
Solid lines represent the model with reservation (car-sharing) while dashed 
lines represent the model without reservation (bike-sharing). 


In E|, for the corresponding bike-sharing model, the 
minimum 2 /(K + 1) for Pf, is reached for pv = K /2 + A/ p. 
We are also interested analytical expressions for this minimum 
for the car-sharing system but it turns out to be tedious again 
because of implicit equations. Therefore, even if we would 
have proven that this minimum is reached for py = 1, an 
asymptotic expression could only be obtained in two cases: 
light traffic (A/p —> 0) or heavy traffic (A/p —> +oo). 













Property 6 (Light traffic): For any py > 0, when X/p 
tends to 0, 


K K 

s = ^ kpy/Pv + °( 1 ) 
k—0 fc=o 


n = i - 


fV 


PV 


pZ +1 ~ i 


+ o(l) 


Proof: The method is standard. It uses that p r = 

-^—TTfrPv~ + o(X/p). The proof is omitted. ■ 

1 Pv I 1 _ __ 

Property [6] means that, in case of light traffic (A /p —> 0), 
the system has the same performance at first order in X/p as 
the homogeneous bike-sharing model . This result coincides 
with the intuition that when the traffic is low, reservation does 
not have time to induce congestion because cars are picked- 
up faster than new users arrive in the system. Particularly, the 
optimal limiting proportion of problematic stations is obtained 
when p v = 1 which gives that Pjf = 2/(K + 1) and s* = 
I \/2 at first order in X/p. This result can be observed on 
Figure [3] where curves for the bike-sharing system get quite 
close to curves for the car-sharing system when the traffic rate 
decreases. 

Property 7 (Optimality for light traffic): If py = 1, as 
X/p —y 0, 


A 


s = 


Pb = 


K K 2 
~2 + 2(K + l) 2 p 
2 2 K X 


K+l (K + l) 3 p 


+ 0((X/p) 2 ), 
0((X/pf). 


Assuming Conjecture |T] the previous quantities are respec¬ 
tively s* and P/. 

Proof: These formulas are straightforwardly obtained, via 
p R = K/(K+l) X/p + 0((X/p) 2 ). U 

Assuming Conjecture Q] Property |7] highlights the differ¬ 
ences of behavior between the bike-sharing system and the 
system with reservation at the optimum. Indeed, s* increases 
less rapidly in the present case and, most of all. If is no more 
constant but has a positive term in X/p when traffic is light. 
However, increasing the capacity I\ dramatically reduces the 
effects of this term. 

Property 8 (Heavy traffic): For all e > 0 and M > e, for 
any py £ [e, M], as X/p —> +oo, 


S = K - J—yJp/X 
V PV 

+ 0(p/ A) 3 / 2 ) 

'A 


3 + py — 2 K 
2 pv 


p/X 


P b = l-(K- 1) 


(py + l)(K-l) 
Vp^k 


+ 0 1 (- 


— 2 > 



Proof: We proof fisrt that pn = y/Kpy^/X/p — (1 + 

pv)/2+ {5 /s - k/ 222^ )+Pv/4 Vp/^+°(p/ x )- The method 
is standard. ■ 


The behaviour of If is interesting in Property [8] as it 
happens to tend to 1 rather quickly compared to the conver¬ 
gence of s and independently of py (the first two terms do 
not depend on py). This is also the reason why we chose to 
dig further in calculations than in Property [6] We shall also 
mention that this speed of convergence is higher because of 
the term p^/(K\Z) which does not exist in P & + . Assuming 
Conjecture [I] to be true, it is interesting to consider the case 
Pv = 1 - 

Property 9 (Optimality for heavy traffic): If py = 1, as 
X/p —> +oo. 


s = K- 

P b = 1 - 



A' - 1 
X/p 


2 — I< 


O 


X/p 
2(K — 1) (X 


s/K 


P 





Assuming Conjecture |T] the previous quantities are respec¬ 
tively s* and If. 

Proof: We obtain these formulas straightforwardly from 
Property [8] ■ 


Problematic stations at a given s, with K = 4. 



Traffic A//./ 


Fig. 4. Evolution of the proportion of problematic stations with traffic (s 
fixed). Solid lines represent this evolution for different values of s. The dashed 
line represents the evolution if for each \//i the optimal parameter s * is 
chosen. Consistently with Figure [3] the lowest curve depends on the traffic. 


C. Maximizing the number of trips 

The operator of the system, either a private company or 
the city council, may choose parameters differently since it 
has to take into account the cost of running the system. In this 
case, the number of successful trips becomes more relevant as 
a metric because the revenue of the operator is more or less 
directly proportional. While it is just an approximated value, 
we assume in the following that T = A(1 — P b ). 

Let us assume p is fixed. We will study T as a function 
of A which is common in transportation. Let us first study 
the maximum number of trips T* = A(1 — If) that can be 
obtained as a function of A. For that, the previous expansions 



















will be very useful. We use for A small the results for light 
traffic of Property [7| which gives 


which means a linear function of the demand for small A. For 
heavy traffic case, from Property [9] as A tends to oo, 


rjn* 


(K - 1 )p + 


2(K — l)/i 3 / 2 

Vk\ 


+ o 



Figure 0 shows the evolution of the number of successful 
trips as a function of X/p. In practice, it is obtained by solving 
numerically equations (J9]» and ifTTb for each value of A//./. We 
observe the linearity when A//i —>• 0. When X/p -A oo, if s is 
hxed, fiy -A 0 which prevents from using Property [8] (py has 
to be bounded). This is consistent with the fact the solid curves 
do not seem to tend to K — 1 whereas the dashed curve does 
(since it corresponds to py = 1 hxed). Another interesting 
observation is that there is no local extremum (T is strictly 
increasing) which is an argument in favour of Conjecture [2] 
Demand A is a function of the price p of a trip. The profit 
is 


n(s,p) = A(p)(l - Pb{s,p)) - C(s), 

where C is the cost function. Maximizing profit leads to 
choose s such that 


dP b 

ds 


1 

Hp) 


C'(s). 


Since C is supposed to be strictly increasing in s and 
Conjecture Q] to hold, the optimal choice in terms of profit for 
the operator is necessarily reached for some so < s*. Thus, 
even under competition pressure, the trade-off between price 
and number of bikes certainly results in choosing s < s*. 


Successful trips at a given s, with K = 4. 



Fig. 5. Number of successful trips T as a function of traffic \/p for different 
values of fixed s (solid lines) and for the optimal value of s (dashed line). 


V. Extension to a Model with Double Reservation 

Car-sharing systems often offer the opportunity to reserve 
online a car and a parking space in the desired station, a while 
before actually picking-up the car. In this section we present 
a model that adresses this particular demand of travel. 

A. Model Description 

At a given station i, new reservations of car are made at 
rate A. At the same time, the user making this reservation also 
wants to reserve a parking space in some station j. If there 
is no car available in i or no available parking space in j, 
the user is rejected. Otherwise, she waits a time exponentially 
distributed with mean 1/v before coming to pick-up her car. 
Then, the journey towards station j takes an exponentially 
distributed time with mean 1/p. The user returns her car at 
station j and leaves the system. 

A three-dimensional state space (reserved cars, reserved 
spaces and free cars) is unsuited to fully describe this process 
as markovian because the duration of reservation of spaces 
is the sum of two exponentially distributed variables (the 
time to pick-up the car and the time to travel). This leads 
to introducing a fourth variable to distinguish between spaces 
reserved by users not yet travelling and users travelling. As we 
will see in the description of transitions, the Markov process 
associated to this model is quite complicated. Let us make an 
approximation, introduce the following process and discuss 
later the relevance of the approximation. 

Let us denote now x = {(j, k, l, m) £ N 4 , j + k+l+m < 
K} and Y^ k l m (t) the proportion of the N stations where, at 
time t: 


II 

j 

Reserved vehicles 

V?(t) = 

k 

Lree vehicles 

R N {t ) = 

l 

Reserved parking spaces by 
users travelling 

II 

Ho' 

2; 

m 

Reserved parking spaces by 
users not yet travelling 

which can be written 


l,m (^) — 

1=1 

l {vp 

N (t)=j-,Vt N (t)=k-,Rf (t): 


As the system is homogeneous, the process ( Y N (t )) = 
( Yyf, j m {t)) (j,k,l,m)ex * s a Markov process. Suppose the pro¬ 
cess is at state (yj,k,l,m)(j,k,l,m)ex- There are three different 
types of transitions: 

• Reservation. At rate A, a user wants to reserve a car in a 
station i. The station is in state ( j, k , l, m) with probabil¬ 
ity Ny : j k i m . Simultaneously, she chooses an arrival sta¬ 
tion i' in state ( j',k',V,m') with probability yj> k',v ,m’■ 
She succeeds if there is at least one free car in i (k > 0) 
and at least one free space in i' (j' + k' + l' + ml < K). 
Therefore, the transition rate is A Nyj t k,i, m yj',k',l',m' if 
k > 0 and j' + k' + V + m! < K. The car reservation 
causes yj t k,i,m to decrease and yj+i,k-i,i,m to increase, 
both by 1/At. The parking reservation causes yj\k',l',m' 
to decrease and yj\k',l',m'+1 to increase, both by 1/N. 








• Cars picked up. Time between reservation and picking 
up is exponentially distributed with parameter v. For 
stations in state (j, k, l, m ), there are jNyj t k,i, m reserved 
cars. Hence, cars are picked up at rate jvyj t k,l,mN. 
But, at the same time, a reserved parking space changes 
of status (from user not yet travelling to travelling) in 
another random station somewhere else, sampled with 
probability proportional to the number of reserved park¬ 
ing spaces in this station. Thus, it occurs in a station in 
state (j', k !, V , to!) with probability 

m 1 Nyji t k',v ,m> _ to! 

where E y [f? r,Ar ] is the average number of reserved 
parking space for users not yet travelling. Finally, the 
transition rate is vj My ^r, N] Nyj,k,l,mVj' , fc 'In the 
original model, this change of status occurs in the very 
same station as the one chosen during reservation: the 
approximation allows to keep the memoryless property. 
The transition causes yj t k,i,m to decrease and yj~i t k,l,m 
to increase. The change of status in the arrival station 
causes yj> t k> ,v ,m' to decrease and yj> t k \to in¬ 
crease. All by 1/N. 

• Cars returned. Considering stations in state ( j,k,l,m ), 
the number of reserved parking spaces is INyj.k.i- As 
trip are exponentially distributed with mean 1/p, cars are 
returned at rate plNyj t k,i■ When a car arrives, it causes 
Vj,k,i,m to decrease and yj,k+i,i-i,m to increase, both by 
1 /N. 

B. Probabilistic Interpretation 

The simplification introduced in the second transition 
allows to write down a differential system of equations similar 
to what is done for the simple-reservation model: the limiting 
empirical distribution y(t) of the stations evolves in time as 
the distribution of some non-homogeneous Markov process 
on x, whose jumps are given by T y (t), updated by the current 
distribution y(t). This process can be seen as the queue length 
vector of four coupled queues with an overall capacity K (see 
Figure [bji. 
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Fig. 6. A typical station as a tandem of 4 queues with overall capacity K. 


As it is done for the simple-reservation model, this 
interpretation leads to 4 equations on the parameters of the 
queues (pr,pv,pr,pv), and a fifth equation is obtained from 
the expression of the average number of vehicles in the system 
s{pR, PViVRiVv) = s - Simple algebra can reduce most of 


them to the equations studied in the former model via a change 
of variables. Indeed, first, p r and py can be omitted as they 
appear to be directly proportional to pR. Remain 3 equations 
on pr and py. Second, if a = A / p. 


Pr= (y 1 + 2 pr 
and a = ^1 + 2—^ a 

leads to pr and py satisfying the exact same equations 
0 and [TOj with parameter a. However, regarding the relation 
with the average number of vehicles, the result is not as simple 
since after the change of variable we have 


s(pv,Pr) 
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E 

k+l<K 



i + p/y \ 

1 + 2ylt /V J 


k Pr 
P V n ■ 


( 22 ) 


The increase of s with respect to its parameters still has 
to be proven because of this new coefficient (k + l ) ■ 

But if we assume 1/V to be sufficiently small, an argument of 
continuity allows to use Theorem [2] Therefore, existence and 
mostly uniqueness of the equilibrium point are straightforward 
with the different results we have for the simple-reservation 
model. 


C. Consequences on Performance 

Thanks to the change of variable, we also have that 
Pb writes the same way as a function of py and pr, with 
parameter d. As a consequence, we still have that a local 
minimum of Pb is reached for py = 1 as mentioned in 
Section |IV] and the value of this minimum is the same as 
for the first model with a traffic ofa = (l + 2^). Similarly, 
Pfc(O) = 1 and Pb —> 1 when py —y +oo. 

However, Equation [22] shows that for some number 
of cars s, the associated solutions pr and py have to be 
higher than the corresponding pr and py in the first model. 
Therefore, if we compare s P&(s) to its analogue in the 
first model with traffic a, it is smaller for low values of s, 
reaches a minimum before and then is higher for s closer to 
K. 


D. Relevancy of the model 

This model introduced above is simplified in the de¬ 
scription of the transition of users picking up their cars. 
However, given the probabilistic interpretation in terms of four 
coupled queues, we can intuitively expect it to have the same 
stationary behavior as the non-simplified process. To confirm 
this, we performed simulations of the original model where 
people actually go to their reserved parking space instead 
of re-sampling one, and derived the steady-state distribution. 
The empirical performance function resulting from this latter 
distribution s i —> Pb(s) seems identical to the one given by the 
model. 



































VI. Conclusion and Future Work 


In this paper we have investigated the influence of reser¬ 
vation on the performance of vehicle-sharing systems. We 
used a stochastic model and then a mean-field limit to study 
the steady-state behaviour of the system and to compare it 
with former homogeneous models without reservation. Though 
reservation makes calculations more tedious, assuming con¬ 
vergence, we still have theoretical results on the existence and 
uniqueness of the limiting stationary state of one station when 
the system gets large. Techniques are standard with a crucial 
use of probabilistic tools in each step. Because results on 
limiting performance, in particular the optimal performance, 
are difficult to obtain, we gave results in the two cases of 
light and heavy traffic. It appears to be the same as for bike¬ 
sharing systems when traffic is light. But, unlike in bike¬ 
sharing systems, saturation occurs when traffic is high or when 
the average number of bikes per station s is close to the 
capacity. We conjectured and partially proved that the optimal 
choice of s is reached for py = 1, like for bikes, which is 
close to K/2 + X/(2p) in the light traffic case, when X/p is 
very small, but no more in general. In heavy traffic, when A/// 
is large, it is of the order of K — f 

To take into account real possibilities in such systems, 
we also studied a network where users can reserve both a 
car and a space before beginning their journey. We studied an 
approximated model (validated by simulations) and we proved 
this new model is mathematically similar to the previous 
model: all is as if the traffic was increased. 

Some conjectures presented in this paper still have to 
be proved. And similarly to what has been done for bike¬ 
sharing systems in £9, further work might focus on different 
algorithms to improve performance or to generalize the model. 
For instance, in terms of the so called passive regulation, we 
investigated the power of two choice in bike-sharing systems 
and it could also have an important impact on car-sharing 
systems. Or, concerning the regulation by the operator, since 
trucks cannot be used with cars we could look at different 
strategies for ambassadors: they move cars one by one from 
full stations to empty ones which would make the model easier 
to define than for trucks that move several bikes together. In 
terms of generalization, for bike-sharing systems (see El), 
results extend to an inhomogeneous model with clusters of 
stations with the same parameters. The same can be done for 
our model. 
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VII. Appendices 
A. Proof of Property\3\ 

Proof: We will prove that, for all ( pv,p R ) G T, 

d s ds 

' -(pv,Pr) > 0 and -—( py,p R ) > 0 


dpR dp v 

by induction on K. For that, let s be denoted by $k ■ 
By a change of indexes, §k can be rewritten 


K 


sk = 


E ( k + l )Pv if E kp k 

k+l<K k=0 


E P k vi f 

k+l<K 


K 

E Pk 

k =0 


where, by definition, 


k fc _j 


i=0 


Define also for ( k,l ) G \ 


n,k = 


EL 

Pk 


Let k > 0 be fixed. We first show that r k ,k-i is an 
increasing function of both p R and py. Indeed 


Pk 

T/c,/c—1 — 

Pk -1 

thus 


pvPk-i + p k R /k\ p k R 1 

-= Pv + 77 -, 

Pk -i k\ p k -1 


dr k , k -i _ _ p R dpk-i 1 
dp v k\ dp v pl_ 1 ' 


(23) 












Pr 


But, by definition of pk , 

dpk- 


1 


dp v 


fc-l 2k—i— 1 fc-2 2k—i—2 

1 = E = Ec + DAnrfrr- 


Uif = 


v- * — 1 )! 

2=1 V ' 2=0 


Pk _ _ 

^ I] r k,K 

/ j Pk k —0 

/c—0 


(k — i — 2)V 

(24) j s strictly increasing in x £ {pv, Pa}, because rk,K = l/rx,k 
And, using that all terms of the sum in the following ' s str | ct ^ decreasing in x.~ 
equation are positive. 


k\pl- x = k\ Y, Pv +V ~ 2 


2 k—u—v 

Pr 


l<.u,v<k 

k—2 


(k — u)\{k — u)! 


2+1 


k\ 


>EpVPk 1 2 E (fc _2-i + j )!( fc-j )! 


that 


and then 


k\ 
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(k-2-i + j)l(k - j)\ (k-i- 2)! 


2+1 


E th _ o _ 


A:! 
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9py 


> 0. 


Moreover, 


4" 1 


dpR 

Using that 


(fc - l)!pfc_i 


1 - 


Pr dpk-i 
kpk-i dpR 


k -1 


dpk-i 2 P k R l 2 n t -s 

= 2^ PvtttH —7ir( fc - 1 - *) = Pfc-2 


dpR 


2=0 


(fc — 1 — *)! 


and kpk-i > PRPk- 2 , because 

k—2 

PR 




{k-2-i)\k 

n k- i-i 

(fc -1-i)! 


2=0 

k—2 


< 

2=0 

<Pfc~l, 


Let us mention that so is constant with x and that 

§k = (1 — uk)sk-i + Kuki 

which yields that 

dsK , T r ~ \dux , M dsK-1 

= (A - SK- i)-5 -1- (1 - u K )- 


. (25) 


dx 


dx 


dx 


For all j, 0 < j < i + 1, while k — 2—i+j < k , it holds 


Since K—Sk- i > 0, < 1 and > 0, by induction 

we can conclude that > 0 for all K > 1. It ends the proof. 


^ (fc-2-* + j)!(fc- j)! (* — * — 2)!’ 

Plugging in equation (1251) and comparing with equation (l24l >. 
it gives 

> o k ^P k ~ 1 
k-Pk-i > Pr dpy ■ 

Therefore, using (l23l >. it allows to conclude that 

drk.k -1 


(26) 


we can conclude that dr g'p R 1 > 0. 

i 

Consequently, if l > k, rgfc = rgi_i is an 

i=k -\-1 

increasing function of py and pr. This gives that uk defined 
by 


























